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ABSTRACT 


Aerial videos of the surf zone at Sand City, Monterey Bay, CA, were aequired 
using an unmanned aerial vehiele. Videos of 26 individual bores were eonverted to still 
images and were georectified, georefereneed, and post proeessed. The size, shape, and 
evolution of the wave generated foam patterns within the surf zone were analyzed. 
The results were tested against two existing hypotheses of foam pattern generation. 
Obliquely Deseending Eddies (ODEs) and self-organization due to bubble rise. Three 
foam regions within the surf zone were reeognized. The largest region, the foam mat, 
eneompasses nearly the entire surf zone and is deseribed as a mat of foam that develops 
obvious foam holes. The areas of the holes grew in area and elongated with time. 
The results were ineonsistent with the two theories of foam hole generation. The fringe 
region is the most seaward foam region and is marked with cireular foam rings that 
become larger in area and more distinct with time. The fringe region data are contrary 
to observations of ODEs but is consistent with the theory of self-organization due to 
bubble rise. The gap region, located between the plunge point and the splash up 
created by the bore collapse, is marked by horizontal foam tubes oriented in the 
cross-shore direction. The foam tubes are likely created in the convergent region 
between two counter-rotating vortices. 
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I. 


INTRODUCTION 


As waves break near the shore, proeesses oeeur that trap air below the ocean 
surface creating bubbles (Thorpe et al. 1999). The size of the bubble determines how fast 
it will rise to the surface, with the larger bubbles rising fastest. Some of the bubbles break 
as soon as they reach the surface; others remain on or just below the surface, creating 
patches of foam. Directly in front of, above, and behind the wave crest, the foam exists as 
a dense mat with few gaps or holes. Thorpe et al. (1999) noted that due to the nature of 
the turbulent processes in the surf zone, the foam left behind as the bore progresses 
toward the beach morphs into coherent patterns, recognizable to the naked eye. 

Waves break near the shore due to kinematic and dynamic instabilities (Thornton 
1979). If the water at the crest of the wave moves faster than the wave itself, the wave 
will collapse and break; this is the kinematic instability. The dynamic instability develops 
when the change in the vertical velocity at the sea surface exceeds gravity. When the 
dynamic instability occurs, the wave becomes too steep (the ratio of the deep-water wave 
height to the wavelength is greater than 1/7), and the wave becomes unstable and breaks. 
The surf similarity parameter, ^b, is widely used as a guide for determining the different 
types of breaking that occur (Battjes 1974). The surf similarity parameter is controlled by 
beach slope angle, breaker height and deep-water wavelength. There are three types of 
breaking waves: spilling (^b<0.4), plunging (0.4<^b<2), and surging (^b>2). Determining 
the breaker type is important because the type defines how the bubbles and foam are 
generated as the wave breaks. 

Longuet-Higgins and Turner (1974), describe how a spilling breaker generates 
foam, using the “entraining plume” model. According to their model, a spilling breaker 
will break gently at the crest of the wave, trapping air. The air-water mixture will stay at 
the surface because the foamy, air-water mixture is less dense than the water surrounding 
it. The foam patch will begin to grow and ride down the front of the spilling breaker due 
to gravity. More air is entrained into the foam patch along the front edge as the patch 
rides down the wave, while more water is entrained into the foam patch from below, 

along the foam-water interface. The foam damps the growth of the crest, and when the 
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damping process exceeds the shoaling process, the crest of the spilling breaker rounds 
and the foam patch slides backward over the crest. Turner and Turner (2011) also discuss 
how foam is generated behind a spilling breaker. Their theory originates with the 
“entraining plume” model, but they add that foam generated in front of the crest can be 
entrained downward and will pass under the crest. Because the foam is less dense and 
more buoyant than the surrounding water, it will float to the surface and will join up with 
the foam that passed over the crest. This foam patch remains behind as the bore continues 
to progress toward the shore. 

As opposed to the gentle breaking process of a spilling breaker, a plunging 
breaker becomes vertical and then overturns, breaking kinematically and throwing a jet of 
water forward of the wave (Peregrine 1983). Derakhti and Kirby (2014) outlined three air 
entrainment mechanisms for plunging breakers. First, as the jet crashes back to the 
surface, the tube closes and entrains air, creating bubbles and foam. Then, as the jet 
crashes into the surface, splashes are sent up forward of the bore. As the bore collapses, a 
cloud of bubbles and water are thrown skyward. As this secondary splash falls back to the 
surface, even more air is entrained in the water, and more bubbles and foam are 
produced. This process of jet formation and splash down can be repeated several times as 
the wave moves toward the shore. Once the wave becomes a bore, more air is entrained 
in a process similar to that described above for a spilling breaker. 

Turbulent processes in the wave-breaking region are not well understood. In the 
past, hot-film and laser-Doppler anemometers were used to measure turbulence in the 
surf zone; however, these point measurements were not sufficient to describe the 3- 
dimensional aspects of turbulence (Ting and Reimnitz 2015). Foam and bubbles in the 
surf zone can cause signal dropout, further complicating efforts to capture information on 
turbulent processes (Nadaoka et al. 1989). In recent years, particle image velocimetry 
(PIV) and volumetric three-component velocimetry (V3V) have been used to measure 
breaking wave generated turbulence in laboratory wave tanks (Ting and Nelson 2011, 
Ting 2013, Ting and Reimnitz 2015). However, these techniques have yet to be adapted 
to field measurements. Nadaoka and Kondoh (1982) observed that the intensity of 
turbulence at the bottom of the water column is correlated with the turbulence at the top 
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of the water column. Therefore, the foam on the surface can be used as a tracer to 
describe total water column turbulence. These turbulent properties are important to 
categorize because they drive processes like sediment transport, water clarity, and the 
transport of pollutants in the surf zone. The presence of foam affects the ability to see 
through the water column as well as the performance of acoustic and optical sensors. 
Understanding how the foam evolves will help the operator better exploit the 
environment. 

Thorpe et al. (1999) and Thorpe (2007) describe two mechanisms for generating 
foam holes in the surf zone; obliquely descending eddies (ODEs)/hairpm vortices, and 
self-organization by rising bubbles. In a wave tank, Nadaoka et al. (1989) observed that 
in the wave breaking region two types of eddies develop, namely horizontal eddies and 
ODEs. The horizontal eddies, which are located around the wave crest, have a horizontal 
axis aligned parallel to the wave crest. The ODEs form behind the wave crest and in the 
region between two horizontal eddies. These ODEs (Eigure 1) have an axis that is aligned 
perpendicular to the wave crest and tilted downward in the offshore direction (Nadaoka et 
al. 1989). Nadaoka et al. (1989) also observed that the ODEs can reach the bottom and 
are a mechanism for suspending sediment in the water column. As the eddy reaches the 
bottom and “kicks up” sediment into the water column the buoyancy of the bubbles in the 
eddy help to further move the sediment upward in the water column. Hairpin vortices, 
also known as reversed horseshoe vortices, have a similar structure to ODEs. The legs of 
the hairpin vortices are made up of two counter-rotating ODEs connected at the base by a 
transverse vortex (Watanabe et al. 2005, Earahani and Dalrymple 2014, Ting and 
Reimnitz 2015). 
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Figure 1. Schematic of Eddies in the Wave Breaking Region 



Horizontal eddies depieted under the breaking wave and alligned parallel to the wave 
crest. Obliquely descending eddies depicted behind the wave, extending from the surface 
to the sea bed. Arrow depicts direction of wave propogation. Source: Nadaoka K., Hino 
M.,Koyano Y., 1989, Structure of the turbulent flow field under breaking waves in the 
surf zone, J. Fluid Mech., 204, 359-387. 


The second possible mechanism for foam hole generation is the process of self¬ 
organization by rising bubbles. Derakhti and Kirby (2014) showed by numerical model 
that in both spilling and plunging breakers the drag force between the bubbles and the 
liquid surrounding the bubbles was the primary mechanism for momentum exchange in 
the vertical direction. Therefore, as the larger bubbles rise rapidly to the surface they 
impart upward momentum on the water surrounding them (Thorpe et al. 1999). When the 
larger bubbles reach the surface, the upward-induced motion continues briefly, causing 
divergence at the surface and a flow of smaller bubbles radially outward, creating a 
region of no foam or bubbles. They also noted that there is a downward motion (due to 
continuity) that occurs at a greater radial distance. 

The primary goal of this thesis is to study the nearly circular patches of foam and 
foam holes that emerge behind bores inside the surf zone. The foam rings that exist 
between the foam mat and the open water will also be studied along with the foam tubes 
that develop between the plunge point and the secondary splash point. The spatial and 
temporal scales of these foam patterns have not been quantified. With the advances made 
in small, inexpensive, unmanned aerial vehicles, persistent overhead imaging of the 


4 





natural surf zone is now possible. Aerial imagery of the surf zone is analyzed to 
determine how the shape and size of the foam holes evolve with time. These foam 
patches are used as a passive tracer to describe the turbulent processes that occur due to 
wave breaking, and the statistics collected from analyzing these foam patches are used to 
quantitatively describe the turbulence. 
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II. METHODOLOGY 


A. EQUIPMENT AND EXPERIMENT DESIGN 

The DJI Phantom 2 Vision+ (hereafter referred to as the Phantom 2) quadcopter, 
unmanned aerial vehiele, UAV (see Figure 2) was used to aequire aerial imagery of the 
surf-zone. The lightweight, quad rotor, battery powered aircraft was chosen for this 
project because of its ease of flight, image capabilities owing to its ability to hover in a 
station-keeping mode, and relatively low cost ($lk). The Phantom 2 has approximately 
20-minute flight time. The Phantom 2 supports a gimbaled, fish-eye lens camera located 
on the underside of the vehicle. The camera can record high definition video at 1080p30, 
1080 horizontal lines of vertical resolution, progressive scan, where every horizontal line 
is drawn in sequence and an image capture rate of 30 frames per second. The camera can 
also continuously record still images at a resolution of 4384x3288 pixels at 1 frame every 
3 seconds. Video mode rather than still image mode was chosen to maximize the 
frequency at which images were collected. The user can choose between either 85° or 
110° recording field of view angles (http://www.dii.com/product/phantom-2-vision- 
plus/specT Based on camera specifications, the field of view and resolution in 
pixel/meter were calculated for various flight altitudes assuming the UAV was directly 
above the surf zone (Figure 3). The FAA maximum flight altitude for unmanned aerial 
systems is 152.4 m (shown as horizontal black line in Figures 3c and 3d) 
lhttp://www.faa.gov/news/press releases/news story.cfm?newsld=l82951 . The field of 
view was decreased by 20% to account for the portion of the image that is lost during 
image processing (image processing is discussed in greater detail below). Based on these 
calculated values, a flight altitude between 76 and 91 meters was chosen to observe the 
estimated width of the surf zone 0(100m), conform to FAA regulations on UAVs and to 
resolve the features of interest. It is desirable to acquire aerial imagery from directly 
above the surf zone with no tilt to the camera to minimize image distortion. However, in 
order to capture the entire width of the surf zone (approximately 100 meters) at a high 
enough resolution and capture all of the ground control points (GCPs), some of which 
were located on the beach, it was necessary to apply a slight tilt to the camera. The 
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wireless connection between the quadcopter and the remote control also limited the 
maximum flight altitude due to signal drop out at altitudes above approximately 100 
meters. The calculations of pixel resolution determined the size of the GCPs required to 
attain sub-pixel accuracy. 


Figure 2. DJI Phantom 2 Vision+ 



The DJI Phantom 2 Vision+ Quadcopter used during this experiment equipped with a 
small 14 megapixel gimbaled camera, compass, live video stream downlink. Source: DJI, 
2015: DJI Phantom 2 Vision +. Accessed 8 October 2015. [Available online at 
http://www.dji.com/product/phantom-2-vision-plus/feature.] 


Figure 3. Resolution (pix/m) and Field of View (m) versus Elevation (m) 
a) b) 






a) resolution versus elevation for a view angle of 110° b) resolution versus elevation for a 
view angle of 85° c) FOV versus elevation for a view angle of 110° (solid line) FAA 
Maximum allowed flight altitude (dashed line) d) FOV versus elevation for a view angle 
of 85° (solid line) FAA Maximum allowed flight altitude (dashed line) 
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Ground control points (GCPs) are needed to georeetify and georeference images. 
For image proeessing purposes it is desirable to have several pixels fall within the GCP, 
this is called sub-pixel aceuraey. Several 0.61-meter diameter plywood cireles were 
painted neon pink and used as GCPs. The pink color was chosen so that the GCPs would 
stand out against the wet sand, water, foam and breaking waves. Five of the GCPs were 
placed atop pipes that were jetted into the sand. One of the pipe GCPs was located within 
the surf zone and broke in half during active wave breaking and had to be removed. 
Three other pipe GCPs were loeated in the swash zone with the last pipe GCP located just 
above the high water line. Four additional GCPs were plaeed on the sand just above the 
high water line and within the field of view (Figure 4). 


Figure 4. Undistorted Image Being Georeetified 



Black boxes indicate search areas, neon green pixels are pixels with color matches to the 
ground control point, black asterisk is ground control point. 
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B. DATA COLLECTION AND PROCESSING 


The video of the surfzone foam patterns was obtained on a mostly dissipative 
beach at Sand City, Monterey Bay, CA on 12 FEB 2015. The beach slope within the surf 
zone is approximately (1:100) (MacMahan et al. 2005). Due to the location of the 
headlands and effects from the Monterey Bay Submarine Canyon, waves at Sand City are 
typically shore normal, resulting in weak alongshore currents and the formation of rip 
currents (MacMahan et al. 2005). Rip currents were avoided in this study. Wave heights, 
wave direction and wave frequencies were determined by analyzing sea surface elevation 
data collected at wave buoy 46042, located just outside of Monterey Bay (Figure 5). The 
wave height on 12 February was 1.87 meters and the peak period was 13.79 seconds. The 
swell wave direction is from the west. The day of collection was specifically chosen due 
to the long wave period and the narrow frequency band of the incoming waves (Figure 5). 
A longer wave period was desirable because it allows more time for the foam and foam 
holes to develop, move and change before the next wave broke. The narrow frequency 
band indicates that the breaking waves are swell waves, and that they are not 
contaminated with waves of other frequencies. In this data set, the observed waves break 
initially as plunging breakers, but quickly become spilling breakers and only one 
jet/splash down cycle is recorded. Eleven videos were collected, with video length 
limited to less than 2 minutes to make the data size more manageable. Continuous in 
flight video downlink was enabled through the use of the DJI Phantom cell phone 
application and the video feed was constantly monitored on a cell phone to ensure that 
the ground control points were in the center of the field of view. The locations of the 
GCPs, with horizontal accuracy on the 0(cm), were determined by post-processing GPS 
data from a survey-grade GPS receiver deployed on the center of the GCPs. 
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Figure 5. Sea Surfaee Elevation Speetrum Buoy 46042 
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Narrow-band sea surface elevation spectrum at buoy 46042 located just outside Monterey 
Bay on 12 February. The wave height was 1.87 meters and the peak period was 13.79 
seconds. The wave direction for the peak period is from the west. 


The video was downloaded and eaeh frame was eonverted to a still image. Every 
10th image was selected, providing 3 frames per second instead of 30. This was done to 
keep the data size more manageable while keeping the image collection frequency high 
enough to capture the features of interest. To remove fish eye lens distortion, the 
camera’s intrinsic parameters must be known or calculated. MATEAB provides a simple 
camera calibration application that determines intrinsic camera parameters using the 
methods described in Heikkila and Silven (1997) and Zhang (2000). These camera 
parameters are then applied to each image to remove the fish eye distortion. 

The undistorted photos were georectified and georeferenced similar to a process 
described by Brouwer et al. (2015). Georectification is the process of taking a two- 
dimensional (2D) image and warping it, using a transformation matrix, to match an earth- 
based three-dimensional (3D) coordinate system. Elsing a projective transformation 
matrix, the image pixel containing the ground control point was matched to the latitude 
and longitude of the ground control point measured in the field. To automate the process, 
a set of search boxes in the image, each containing a ground control point, were created 
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(outlined in black in Figure 4). It is necessary to use a search box because the UAV 
moves slightly with wind gusts. Each pixel within each search box was then evaluated to 
determine if the pixel’s red-green-blue color value fell within a range of values matching 
that of the pink ground control point (all matching pixels displayed in green in Figure 4). 
The center of mass of each set of matching pixels was calculated and that pixel was then 
stored as the image pixel location corresponding to a ground control point (shown as a 
black asterisk in Figure 4). The process of georectification was essential because 
knowing the size of each pixel was required in order to determine the size of the foam 
and foam holes. Georeferencing allows the user to spatially compare georectified images 
between frames. Each pixel in an image is assigned a corresponding x and y reference 
value, thus enabling the conversion from pixel space to latitude and longitude space. 
Once the images were undistorted, georectified and georeferenced, they were rotated so 
that the bore would be parallel to the x-axis (Figure 6). This was done so that the analysis 
is described in cross-shore and along-shore directions relative to the bore. 
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Figure 6. Georeferenced, Georectified and Rotated Image 
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Rectangle A depicts the fringe region, rectange B depicts the gap region, rectangle C is 
the plunging breaker, circle D represents a foam hole, rectangle E depict a top box, 
rectangle F depicts a bottom box, rectangle G depicts a self similar bore, W depicts the 
width of the surfzone. The yellow circles indicate the GCPs and the pink circle shows the 
UAV operator. 


The undistorted, reetified, refereneed and rotated images were then individually 
analyzed and only a small portion of eaeh image, the part filled with foam and holes, was 
saved for further analyses. As a wave broke, the area just behind the wave was seleeted as 
the first top box (Figure 7a). The left and right bounds of the box were seleeted so that 
any wave refleetions from the shoreline or rip eurrents in the image would not be 
ineluded in the sub-box. As the bore progressed toward the shore in eaeh subsequent 
image, the seleeted box beeame larger in the eross-shore, extending from either the edge 
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of the foam or the front of the next bore on the seaward side to the baek edge of the bore 
on the shoreward side. 

As the bore progressed toward the shore there was a point where the foam holes 
developed two distinet regimes (Figure 6 box E and F and Figure 7b). The more seaward 
holes had grown and were now signifieantly larger than the “new” more shoreward holes. 
At this point the box was split into two boxes depending on the size of the holes. The top 
box is loeated between either the front edge of the next bore and the top of the bottom 
box (Figure 6 box E) or between the seaward side of the foam and the top of the bottom 
box (Figure 7b box A). The bottom box was loeated further shoreward, between the 
bottom of the top box and the baekside of the bore (Figure 6 box F and Figure 7b box B). 
The top and bottom boxes were allowed to grow and change as the bore progressed 
toward the shore and represent different depths. As a new wave broke the process of 
analyzing the image and selecting sub-boxes started all over again. Twenty-six individual 
bores and their associated top and bottom sub-boxes were stored for further analysis. 

Figure 7. Selection of Top and Bottom Sub Boxes 


a) b) 



0 50 100 150 0 50 100 150 

local X [m] local x [m] 

Georectified, georeferenced and rotated image depicting a top and bottom sub-box 
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These subset images were then eonverted to gray seale images, pixel value 
between 0 and 1, and then into binary images, pixel value of either 0 or 1. To eonvert the 
images into binary images a threshold value had to be determined a priori. Any pixel that 
had an intensity value larger than the threshold value was eonverted to a 1 (white) and 
any pixel that had an intensity value less than the threshold value was eonverted to a 0 
(blaek). Initially a 0.5 threshold value was applied to all the images. However, after elose 
inspeetion, it was found that using only one threshold value for all the images did not 
appropriately elassify the holes and foam, and some of the important data were lost. To 
aecount for ehanges in image intensity over a set of images, Carini et al. (2015) and 
Podobna et al. (2010) ereated histograms of image intensity to determine the appropriate 
threshold value to use for eaeh image. In a manner similar to this approach, the mean 
intensity value of each image was calculated and used as the threshold value. Again, this 
mean value method was not sufficient to appropriately represent foam and hole. 
Therefore, a manual approach was adopted. 

To manually determine the threshold value, an interactive selection technique 
(Grande 2000) was performed. For each second (every third image), a subset gray-scale 
image (Figure 8a) from the longest dataset of top boxes was displayed on the screen 
along with five binary images, each with a different threshold value (Figure 8). In this 
example, the image displayed is 13 seconds after wave breaking, and the threshold values 
displayed on the image are; 0.475, 0.45, 0.425, 0.4 and 0.375 (Figure 8). If the threshold 
value is set too high, gray pixels that should be turned white are turned black, which 
causes several holes to become connected and appear much larger than they actually are. 
If the threshold value is set too low, the gray pixels that should connect to other gray 
pixel and form a hole turn white, making the holes appear much smaller than they 
actually are. The threshold value that best represented the foam and holes was recorded. 
In this case the threshold value selected was 0.475 (Figure 8b). When viewed on a 
computer screen, it is much easier to differentiate between the different threshold values 
and their impact on the image. A third-degree polynomial was fit to the threshold values 
for all times using the longest bore record (Figure 9 grey line) and was then applied to all 
the bores. The threshold determination was similarly estimated on the longest dataset of 
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the bottom boxes and another third-degree polynomial equation of threshold values was 
reeorded (Figure 9 blaek dashed line). Note that the threshold values were different 
between the bottom and top boxes. These threshold values were then used to convert the 
gray-scale images into binary images for further analysis. The black regions in the binary 
image correspond to the foam holes and the white regions correspond to the foam. 


Figure 8. Determining the Threshold Value 



local X [m] local x [m] 


a) gray-scale subset image b) binary image with threshold=0.475 c) binary image with 
threshold=0.45 d) binary image with threshold =0.425 e) binary image with threshold 
=0.4 f) binary image with threshold =0.375 
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Figure 9. Third-Degree Polynomial of Threshold Values 

for Top and Bottom Boxes 



Time series of threshold values for determining whether a pixel should be a 0 or a 1 
based on the interactive selection technique and a third-degree polynomial best-fit line. 
Gray solid line is used for the top boxes and the black dashed line is used for the bottom 
boxes. 


The boundaries between the blaek regions and white regions were determined and 
used to ealeulate the area of the foam holes. While eondueting a quality eheek of the data, 
it beeame apparent that some of the regions that were identified as foam holes were in 
faet small unimportant features within the image, artifaets of the eonversion to a binary 
image. A minimum area was determined that removed the unimportant features while 
maintaining the features of interest. To determine this minimum area, four of the same 
binary images were displayed with the boundary between the blaek and white pixels 
displayed in red and green (Figure 10). The red boundaries show the areas that are 
smaller than the assigned eutoff area and the green boundaries show the areas that are 
larger than the assigned eutoff area. By using this teehnique, a eutoff area of 0.35m was 
ehosen and any blaek region with an area less than this eutoff was disregarded and 
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removed from further analysis. The foam holes that are not outlined by anything are 
connected to the border and are also removed from further analysis. 


Figure 10. Determining the Minimum Allowable Area 





20 30 

local X [m] 


Features included in green, features ignored in red a) cutoff area = 0.4m^ b) cutoff area = 
0.35m^ c) cutoff area = 0.3m^ d) cutoff area = 0.25m^. Features not surrounded by a color 
are connected to the border and are ignored. 


To determine how the foam holes changed with time, various shape factors were 
analyzed. Shape factors are non-dimensional quantities calculated from the dimensions of 
the shape, which describe the shape of an object relative to another standard object 
(Wojnar 2000). The shape factors used in this analysis are elongation, circularity, and 
compactness. Elongation, defined as. 


Elongation - 



( 1 ) 


V w 

is a measure of the largest length of the shape, 1, versus the shortest width of the shape, w. 
An elongation value of 1 represents a perfect circle and an elongation value much less 
than one represents a stretched ellipse. Both circularity and compactness values were 
evaluated but did not show appreciable change and will not be discussed further. 


18 









A similar process to the one deseribed above, was applied to the fringe region 
between the foam mat and the open water to determine the shape, size, and evolution of 
the foam patterns in the fringe region (Figure 6 box A). This method of analyzing the 
foam patterns in the fringe region was unsuceessful beeause the primary foam rings were 
filled with less dense foam that impaeted the ability to determine a threshold value that 
appropriately classified foam from no foam regions. In order to determine the shape and 
size of these foam patterns, the cross-shore and along-shore diameter of the foam rings 
were caleulated manually, and from these diameters, the areas and elongations of the 
foam rings were estimated. 

The largest wave breaking events produced a small gap between the plunge point 
and the splash up due to bore collapse (Figure 6 box B). In this gap, small horizontal 
foam tubes developed in the cross-shore direetion, eonneeting the plunge point and the 
location of the splash up due to bore collapse. To charaeterize the size and separation of 
these foam tubes, the gap region was seleeted and eonverted to a binary image. A line 
was selected in the along-shore direction that bisected these foam tubes. This biseeting 
line was then used to determine the width and spaeing of the foam tubes as a function of 
time. 
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III. RESULTS 


An initial visual analysis of the video and still images was eondueted to determine 
the general pattern of foam and foam hole distribution as the plunging wave broke and 
progressed toward the shore as a self-similar bore (Figure 6). There appears to be three 
distinct foam pattern regions within the wave breaking/run-up cycle. The first of these 
regions, henceforth called the gap region, occurs only during the largest plunging waves 
and is located between the plunge point and the location of the splash up due to bore 
collapse behind the wave front (Figure 6 box B). In this gap region there are several 
horizontal foam streaks oriented in the cross-shore direction that seem to connect the 
primary bore to the secondary splash point. This gap region exists for 1-4 seconds, 
depending on the height of the wave and the distance between the plunge point and the 
secondary splash. The next region of interest, hereafter called the fringe region, exists 
seaward of the wave immediately after the wave breaks (Figure 6 box A). A consistent 
line of boil like, foam rings, aligned in the along-shore direction mark the fringe region. 
These circular boils appear to grow slightly and become much more pronounced with 
elapsed time from wave breaking. The boils within the fringe region appear more 
elongated if the wave breaks closer to the shore or in a region still covered by the 
previous waves’ foam mat. The third region is the primary foam region and it can be 
described as a dense mat of foam that appears directly behind the self similar bore just a 
few seconds after the wave breaks (Figure 6 boxes E and F). This mat of foam grows in 
extent as the bore progresses toward the shore. Several more seconds after the wave 
breaks the dense mat of foam begins to develop a visible pattern of nearly circular foam 
holes (Figure 6 circle D) starting first at the most seaward extent of the foam mat and 
extending toward the shore as the wave propagates. These foam holes appear to grow and 
become less circular (more elongated) with time elapsed since wave breaking and 
distance from wave crest. 
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Twenty-six bores are evaluated to quantitatively deseribe the evolution of the 
foam holes. The areas of the foam holes are ealculated for eaeh bore at eaeh time step for 
the top sub-boxes. A mean area is then ealeulated for each of the twenty-six bores at each 
time step. These mean areas are then averaged to provide a time series of the ensemble 
mean area of all the bores (Figure 11a). The 95% confidence intervals and the linear 
regressions are also plotted. Each bore exists for a different length of time. Therefore, 
near the end of the time series there are fewer data points to average, causing the 
confidence interval to spread. The data are truncated past 16.33 seconds to account for 
this spread and the linear regression is only plotted for the first 10 seconds. This same 
process is repeated for the bottom sub-boxes (Figure 11b). For the top boxes the mean 
area of the foam holes is initially between 1 and 1.4 meters . For the first 4 seconds the 
area of the foam holes decreases linearly at a rate of 0.08 meters /second with a minimum 
mean area of approximately 0.9 meters reached at t=4 seconds. After t=4 seconds the 
area of the foam holes increases linearly at a rate of 0.11 meters /second, reaching a 
maximum mean area of 2.3 meters at t=15.6 seconds. For the bottom boxes the mean 
area of the foam holes at t=0 seconds is between 0.7 and 0.8 meters . The mean area 
increases linearly at a rate of 0.084 meters /second and reaches a maximum value of 1.86 
meters at t=9.3 seconds. The bottom boxes do not show the initial decrease in area that 
was seen in the top boxes because by the time the wave is in the bottom box region it is 
not a plunging breaker but rather a self-similar bore and there is no splash down cycle to 
cause the large foam holes. 
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Figure 11. Area of Foam Holes All Bores (m ) versus Time (see) 




a) Top boxes. Slopes of best fit lines: -0.081 with an value of 0.78 and 0.11 with an 
value of 0.94. b) Bottom boxes. Slope of best fit line 0.084 with an R^ value of 0.81. Gray 
line is 95% eonfidenee interval, blaek dashed line is the best fit line, blaek vertieal line is 
the peak wave period. 

A similar proeess as was eondueted above to ealeulate mean areas is used to 
caleulate the mean elongation (Equation 1). A time series of these mean values is plotted 
to determine how the elongation of the foam holes ehanges with time (Figure 12). The 
elongation of the foam holes for the top boxes (Figure 12a) is initially 0.94 and deereases 
exponentially at a rate of 0.018/seo. The minimum elongation value at the end of the time 
series (exeluding the final point) is 0.73 at t=16 seeonds. The elongation of the foam 
holes for the bottom boxes at t=0 is 0.91 and deeays exponentially at a rate of 0.021/seo. 
The minimum mean elongation of the foam holes in the bottom box is between 0.76 and 
0.77 at t=9.6 seeonds. 
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Figure 12. Elongation of Foam Holes versus Time (sec) 



a) Top boxes. Exponential deeay of best fit line -0.018/see with an value of 0.83. b) 
Bottom boxes. Exponential deeay of best fit line -0.021/see with an value 0.72. Gray 
line is 95% eonfidenee interval, blaek dashed line is the best fit line, blaek vertieal line is 
the peak wave period. 


While the time series of mean areas and mean elongations provide a general 
understanding of the foam hole evolution, they are not sufficient to describe the small 
details of the evolution of individual foam holes. A visual analysis of a small section of 
both top and bottom sub-boxes yields some interesting results that describe even more 
detail on how the foam holes evolve (Figures 13, 14 and 15). The areas of some of the 
foam holes grow by merging with other foam holes (Figure 15 red circles and pink ovals. 
Figure 13 pink ovals) while the area of certain other foam holes grow independent of 
their surroundings. Foam circles also appear to split apart (Figure 15 blue squares), thus 
decreasing in size. The elongation statistics discussed above show that the foam holes 
stretch but do not provide information on their axis of dilatation. By analyzing the small 
sub-boxes, it is apparent that the primary axis of dilatation is in the cross-shore direction 
(Figure 15 pink ovals. Figure 13 red circles and blue boxes. Figure 14 red circles). The 
statistics above do not describe how the foam holes move in reference to a fixed 
reference frame. By analyzing the small sub-boxes, it appears that the foam holes 
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translate slowly toward the shore (Figure 13 red circles and blue boxes, Figure 14 red 
circles). The foam holes appear to migrate faster in the bottom sub-boxes than they do in 
the top sub-boxes. There does not appear to be any rotation of the foam holes over time. 
The elongation of all the foam holes for one bore was plotted against the area of all the 
foam holes for the same bore. The elongation and area are not statistically correlated. 


Figure 13. Zoomed in Region of Bottom Boxes 



local X [m] local x (ml 


Time series of zoomed in region of bottom boxes for one bore starting at a) t = 4.33 
seeonds with time interval of .33 seeonds. Red eireled region depiets the elongation of a 
foam hole in the eross shore direetion and the transport of the foam hole toward the 
shore. Blue squared regions show elongation and growth of foam holes in the eross shore 
direetion as well as transport in the shoreward direetion. Pink eireled regions depiet 
merging of foam holes and a small pateh of remnant foam within a foam hole. 
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Figure 14. Zoomed in Region of Bottom Boxes 



Time series of zoomed in region of bottom boxes for one bore starting at a) t = 6.66 
seconds with time interval of .33 seconds. Red circled region depicts the elongation of a 
foam hole in the cross shore direction, the transport of the foam hole toward the shore 
and the growth of the foam hole. 
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Figure 15. Zoomed in Region of Top Boxes 
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Time series of zoomed in region of top boxes for one bore starting at a) t = 10.66 seeonds 
with time interval of .33 seeonds. Red eireled region depiets the growth of a foam hole by 
merging with other foam holes. Blue squared regions show foam holes splitting apart. 
Pink oval regions depiet merging of foam holes and elongation in the eross-shore 
direction. 


27 







There were eight breaking waves that produeed pronouneed fringe regions (Figure 
6 box A). The areas of five foam rings along the fringe region were ealculated at each 
time step for each of these eight breaking waves. The mean area of the foam rings at each 
time step for each breaking wave was computed and then used to calculate a time series 
of the mean ensemble area of the foam rings for all eight of the breaking waves (Figure 
16). The 95% confidence interval and linear regression were also calculated and are 
plotted along with the mean ensemble area. The mean area of the foam rings is initially 
5 meters and increases at a rate of 1.05 meters /second. After about 6 seconds the foam 
rings stop growing, maintaining a constant area of about 12 meters^. Initially, the foam 
rings within the fringe region are filled with tiny bubbles. As time elapses, the center 
of the foam rings appear dark black, indicating there are no more bubbles present. The 
elongation of the foam rings was calculated but did not exhibit appreciable change. 

Figure 16. Area of Foam Rings in Fringe Region (m ) versus Time (sec) 



Slope of best fit line 1.05 meters^/seeond with an value of .96. Gray line is 95% 
eonfidenee interval, blaek dashed line is the best fit line. 
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The gap region between the plunge point and the seeondary splash point is 
marked by horizontal foam tubes, oriented in the eross shore direction that seem to 
connect the plunge point to the secondary splash point (Figure 6 box B). Only one wave 
produced the conditions necessary to visualize the foam tubes for sufficient analysis. The 
width of these foam tubes was calculated and a time series of the average width of all the 
foam tubes at each time step was plotted. The foam tubes are only visible from the 
surface for approximately 3.5 seconds and they decrease in width with time. The initial 
width of the foam tubes is approximately 0.25 meters and decreases to a width of 
0.1 meters. 
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IV. DISCUSSION 


Thorpe et al. (1999) and Thorpe (2007) deseribe two physieal proeesses that 
generate foam holes in the surf zone. The first of these proeesses are the ODEs/hairpin 
vortiees that develop as a by-produet of breaking waves. ODEs and hairpin vortiees are 
observed in wave tanks as well as numerieal modeling studies. Observing this 
phenomenon in the natural surf zone is diffieult and to our knowledge no observations or 
measurements of ODEs/hairpin vortiees have been made in the natural surf zone. The 
notion of using surfaee foam as a passive traeer to identify these features, and quantify 
their size and evolution is dependent on the ODE/hairpin vortex penetrating the oeean 
surfaee. Reeent studies on ODEs/hairpin vortiees have been foeused on the layer between 
the free surfaee and above the bottom boundary layer (Ting 2006, Earahani and 
Dalrymple 2014, Ting and Reimnitz 2015) and do not eonsider the effeets of the free 
surfaee on the ODE/hairpin vortex or the effeet of the ODE/hairpin vortex on the free 
surfaee. Ting and Reimnitz (2015) reason that the ODEs/hairpin vortiees must reaeh the 
free surfaee but what effeet these features have on the transfer of mass at the surfaee is 
unknown. 

Several studies on ODEs and hairpin vortiees have been eondueted using 
entrained bubbles as a flow visualization tool to deseribe and traek the evolution of an 
ODE. These studies have been eondueted in wave tanks, and the imagery has been 
aequired from the side to visualize the ODE as it extends down toward the bottom. 
Nadaoka et al. (1989) displayed a plan view photograph from one of their tank 
experiments and they assert that the “spotty features” in the image eould be due to ODEs, 
but there are no data direetly tying the surfaee foam spots to the presenee of an ODE. 
Watanabe et al. (2005) deseribe the proeess of air entrainment in an ODE. They observed 
that bubbles with a eertain amount of vortieity will be eaptured by the ODE and that these 
trapped bubbles will eoneentrate along the axis of rotation of the ODE, where the 
pressure is smallest. If the ODE does penetrate the surfaee, as Ting and Reimnitz (2015) 
elaim, then the eireular pattern on the surfaee, depleting the ODE, would be filled with 
bubbles and not free of bubbles as is evident in the eolleeted imagery. 
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The measured area of the foam holes and the growth rate of the foam holes are 
ineompatible with the hypothesis that ODEs/hairpin vortiees ereate the discemable foam 
pattern on the surfaee is. The areas of the foam holes are between 0.8 and 2.0 meters and 
the diameters are approximately 0.25 to 0.64 meters. The areas are also eonsistent 
between the top boxes, which are located farther offshore, and the bottom boxes, which 
are located closer to the beach. The size of ODEs/hairpin vortices is dependent on the 
depth of the water, where the width of the vortex is approximately 14 of the water depth 
(Ting and Reimnitz 2015). Video data herein shows that the area of the foam holes is 
independent of water depth, contrary to the concept of ODEs/hairpin vortices. Watanabe 
et al. (2015) observed through numerical modeling that the ODE elongates, and therefore, 
stretches over time owing to the conservation of angular momentum, the angular velocity 
increases and the radius of the ODE decreases. As described in chapter three of this 
paper, the foam holes consistently grow with time; therefore the foam holes within the 
dense mat of foam are not a passive tracer of an ODE/hairpin vortex in the water column 
below. 

The second proposed method for foam hole generation is the theory of self¬ 
organization by rising bubbles. Bubbles of various sizes, and therefore buoyancies, are 
entrained into the water column during wave breaking. The larger bubbles rise faster than 
the smaller bubbles, imparting upward momentum on the surrounding liquid as they rise. 
When the large bubbles reach the surface they break and the momentum they imparted on 
their way to the surface spreads out laterally from the center of where the bubble burst. 
The smaller bubbles rise more slowly due to their lower buoyancy value. These small 
bubbles do not burst upon reaching the surface due to their size and surface tension. 
These bubbles are then spread out laterally from the center of their rise location by the 
momentum imparted by the large bubbles. 

There are two pieces of evidence that do not support the theory of self¬ 
organization by rising bubbles. The first is the upward and outward momentum imparted 
by the large bubbles as they rise. Shortly after the large bubbles reach the surface the 
upward and then outward momentum will cease. The video evidence and statistical 
analysis of the foam holes shows that the foam holes continue to grow indefinitely and 
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not just for a few seconds as would be expected from the self-organization argument. The 
second piece of evidence that is incompatible with the self-organization by rising bubbles 
theory is the blackness of the foam holes immediately after formation. According to the 
self-organization by rising bubbles theory, there would be a stream of small bubbles 
rising after the large bubbles that would make the center of the foam hole appear gray or 
white for some time. 

The theory of self-organization by rising bubbles could, however, describe the 
foam rings that develop in the fringe region (Figure 6 box A). These foam rings initially 
appear as bright rings of white foam with lighter, grayer foam in the center. This is 
consistent with the slow rise rate of the small bubbles into the center of the void that was 
produced by the bursting of the large bubbles. The center of the foam rings also become 
blacker and less gray with time elapsed since wave breaking, owing to the outward 
movement of the small bubbles. The rings grow for approximately 6 seconds and then 
remain a constant size, which corresponds with the outward momentum ceasing shortly 
after the large bubbles reach the surface. 

The foam tubes that develop in the gap region of the largest plunging breakers are 
similar to the aerated vortex filaments that were predicted in the numerical simulation 
conducted by Lubin and Glockner (2015) and the surface scars simulated by Saruwatari 
et al. (2009). Much like the gap region foam tubes described in chapter three, the aerated 
vortex filaments modeled by Lubin and Glockner (2015) form between the plunge point 
and the splash-up location and are oriented in a cross-shore direction. These cross-shore, 
horizontal vortices develop in counter rotating pairs. In the upwelling region between a 
pair of vortices, the surface experiences divergent flow and in the downwelling region of 
a pair of vortices, the surface experiences convergent flow (Saruwatari et al. 2009). The 
surface foam formed by the breaking wave is concentrated along the region of 
convergence and downwelling. Because the foam bubbles are buoyant they remain at the 
surface and form distinct foam tubes in the cross-shore direction. In the upwelling region, 
bubble free water is brought to the surface and forms the bubble-free gaps between the 
foam tubes. 
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From the data collected in this experiment, there is still uncertainty as to the 
physical mechanisms that create and evolve the foam patterns in the surf zone. The next 
step in developing an understanding of these foam patterns is to conduct feature tracking 
on the images. In order to feature track the foam patterns, a stationary viewing platform is 
necessary. Feature tracking of the imagery discussed above was attempted, and although 
the georectification and georeferencing processes remove most of the wobble from the 
images, some remains, making the feature tracking process not possible. Another aspect 
that makes feature tracking of the foam patterns difficult is the changing nature of the 
foam. Unlike feature tracking a car or a face, the foam patterns are constantly changing 
and evolving; therefore, high-resolution stationary imagery is necessary. If feature 
tracking is successfully conducted, the surface macro-scale turbulence within the surf 
zone, which organizes the transport of foam at the surface, can be analyzed from a 
Lagrangian point of view. Lagrangian coherent structures (LCSs) are material lines that 
are advected by the flow (Peacock and Haller 2013). LCSs are classified as either 
attractors or repellers. In the region of these attractors and repellers, material develops 
into fillaments, similar to the foam patterns in the surf zone (Figure 17). 

Figure 17. Attracting and Repelling Lagrangian Coherent Structures 


/ = / 




a) an attracting LCS b) a repelling LCS. Source: Peacock, T., and G. Haller, 2013: 
Lagrangian coherent structures: The hidden skeleton of fluid flows. Phys. Today, 66, 


41-47. 
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V. CONCLUSION 


Unique surf zone imagery, aequired from a UAV at Sand City, Monterey Bay, 
CA, is analyzed to deseribe the eoherent foam patterns that evolve within the surf zone. 
On the day of eollection, the narrow band waves were categorized as plunging, shore 
normal, breakers with a period of 13 seconds. Three distinct foam pattern regions are 
observed within the imagery (Figure 6). There is a dense mat of foam that develops into a 
coherent pattern of foam holes that grow and elongate with time. The fringe region exists 
on the seaward edge of the foam mat and is marked by rings of foam aligned in the along 
shore direction. These foam rings become more defined and larger with time. There is 
also a gap region that exists only during the largest plunging waves. The gap region 
forms between the plunge point and the secondary splash and is marked by foam tubes, 
aligned in the cross shore direction that appear to connect the plunge point to the 
secondary splash. 

There are two hypotheses that describe the presence of the foam holes within the 
surf zone. The first hypothesis is that the ODEs/hairpin vortices created by the breaking 
waves cause the spotty foam appearance on the surface (Thorpe et al. 1999 and Thorpe 
2007). The second hypothesis is that self-organization by rising bubbles causes foam 
holes within the surf zone (Thorpe et al. 1999 and Thorpe 2007). By analyzing the UAV 
imagery of the surf zone, both of these theories are tested. 

The first step to analyzing the video imagery is to convert the video to still 
images. Then, the still images are georectified and georeferenced so that the size of the 
features within the surf zone can be measured. The color images are then converted to 
binary images, using an interactive selection technique (Grande 2000), so that the foam 
holes can be distinguished from the foam. The foam mat region is split into two sub¬ 
boxes based on the size of the foam holes (Figure 7). Twenty-six individual bores are 
analyzed and the area and elongation statistics of the foam holes in each box are 
calculated. The mean area and elongation are calculated for each bore and each box, and 
then a time series of the ensemble mean area and elongation is plotted for all the bores, 

for each box. Also plotted are the linear regression and the 95% confidence interval. 
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The mean area of the foam holes for the top box deerease for the first 4 seconds 

and then increases at a rate of 0.11 meters /second. The mean area of these foam holes 

2 2 

range from 0.9 meters to 2.3 meters . The mean area of the foam holes for the bottom 
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box range from .7 meters to 1.86 meters . The mean area of the foam holes for the 
bottom box increases linearly at a rate of 0.084 meters /second. The mean elongation of 
the foam holes for both the top and bottom box was found to decrease, becoming more 
elongated with time. For both the top and bottom box, the elongation values were initially 
1.0 and decreased linearly at a rate of 0.014 for the top box, and 0.017 for the bottom 
box. The final elongation value for both the top and bottom box were between 0.7 and 
0.75. Based on these statistical results and visual analysis of the video data, both 
hypotheses for foam hole generation within the surf zone are incompatible with results. 
The next step in quantitatively describing the foam patterns within the surf zone is to 
conduct feature tracking, which requires a stationary viewing platform. Once feature 
tracking is conducted, the surf zone foam patterns can be viewed from a Lagrangian 
perspective. The LCSs can be calculated to provide information on the material lines that 
concentrate and advect flow, leading to an understanding of the macro-scale turbulence 
within the surf zone. 

The mean area and elongation of the foam rings along the fringe region were 
analyzed, and the rings grow at a rate of 1.05 meters /second. The foam rings do not 
elongate with time. From a visual analysis, the foam rings are initially filled with small 
bubbles, but as time progresses the center of the foam rings appear black indicating all 
bubbles have disappeared. It is reasonable to assume that the generation of foam rings 
within the fringe region is caused by the self-organization due to bubble rise and not the 
presence of ODEs/hairpin vortices. 

Within the gap region, the width of the foam tubes, and the spacing between the 
foam tubes were analyzed. Only one wave produced conditions conducive to foam tube 
generation. These foam tubes are only visible from the surface for a few seconds and 
decrease in width over time. These foam tubes are not likely caused by the ODE/hairpin 
vortex theory or self-organization due to bubble rise, but are rather created by the 
convergent regions between a pair of counter rotating vortices, while the gaps between 
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the foam tubes likely indieate the divergent regions between a pair of eounter rotating 
vortices (Saruwatari et al. 2009 and Lubin and Glockner 2015). 
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